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We consider the contributions to g M — 2 from fourth generation heavy neutral and charged leptons, 
N and E, at the one-loop level. Diagrammatically, there are two types of contributions: boson- 
boson-7V, and _B-_E-boson in the loop diagram. In general, the effect from N is suppressed by 
off-diagonal lepton mixing matrix elements. For E, we consider flavor changing neutral couplings 
arising from various New Physics models, which are stringently constrained by /j, — > &y. We assess 
how the existence of a fourth generation would affect these New Physics models. 
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PACS numbers: 14.60.Pq, 14.60.St 



I. INTRODUCTION 

The fourth generation has been viewed as out of fa- 
vor [1] since a long time, because of electroweak precision 
tests (EWPrT), and neutrino counting on the Z peak. 
However, the severeness of the S parameter constraint 
from EWPrT has been questioned recently [2], while we 
know that the neutrino sector is much richer than origi- 
nally thought because of neutrino oscillations. With the 
advent of the LHC, we now have a machine which can 
discover or rule out the 4th generation by direct search, 
once and for all [3] . Currently, the Tcvatron has set strin- 
gent limits [4] on t' via t' — > qW search. 

It was recently pointed out [5] that the existence of 
a 4th generation could have implications for the baryon 
asymmetry of the Universe (BAU). By shifting the Jarl- 
skog invariant [6] for CP violation (CPV) of the 3 gener- 
ation Standard Model (SM3) by one generation, i.e. from 
1-2-3 to 2-3-4 quarks, one gains by more than 10 13 in ef- 
fective CPV, and may be sufficient for BAU ! Recent de- 
velopments in CPV studies at the B factories [7] and the 
Tevatron [8] suggest the 4th generation could be behind 
some hints for New Physics in b — > s transitions. From 
a different perspective, whether from effective 4-fcrmion 
interactions [9] , or from holographic extra dimension con- 
siderations [10] (the two are complementary), there are 
also recent interest in very heavy 4th generation quarks, 
where their heaviness could be responsible for inducing 
electroweak symmetry breaking itself. 

With renewed interest in the existence of a sequential 
4th generation, and with experimental discovery or refu- 
tation expected at the LHC in due time, we turn to the 
lepton sector. Our goal is modest: if a 4th generation 
exists, what are the implications for the most prominent 
probes with charged leptons, i.e. muon g — 2, /j — * &y, 
and t — > £j? 

The difference between the experimental value and the 
SM3 prediction of muon g — 2 has been around for some 
time now [111. That is, 



bounds on lepton flavor violating (LFV) rare decays, such 
as [11 



B{n -> ei) < 1.2 x 10" 



and the r decay counterpart [12] 

ei) < 1.1 x 10 -7 , 
irf) < 4.5 x 10" 8 , 
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at 90% C.L. These limits could be improved further in the 
near future. The MEG experiment, a ji — ► e^y search ex- 
periment aiming at a sensitivity of 10 -13 [13], has started 
its physics run in 2008. The short 2008 run alone is ex- 
pected to bring the limit below 10 -12 . Though the limits 
on t — ► i/^f from the B factories, Eqs. (3) and (4), will 
soon be limited by B factory statistics, a Super B Fac- 
tory upgrade could push down to the 10 -8 region, which 
become background limited (for outlook, see Rcf. [14]). 

Can the effect of 4th generation leptons show up in 
the probes of Eqs. (l)-(4)? How would these processes 
constrain New Physics models in the presence of a 4th 
generation? In this paper we start our discussion from 
a diagrammatic point of view, and in so doing, correct 
some errors in the literature. 

The 4th generation neutral lepton N can enter the 
loop with charged vector boson or scalar boson H ± , 
which we plot for fifij coupling in Fig. 1(a) for illustra- 
tion. The diagrams for fi — > and r — > £"/ are quite 
similar. The W + W~N loop is controlled by the lepton 
mixing matrix elements, while the H + H~ N loop may be- 
come important because of m^v- The charged lepton E 
can enter the loop with neutral scalar and pseudo scalar 
bosons h and ^4°, or a neutral vector boson Z' , as illus- 
trated in Fig. 1(b). However, these neutral bosons would 
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where a M = (g^ — 2)/2. The difference of over 3.4er has 
aroused a lot of interest. We also have very stringent 



Fig. 1: (a) Boson-boson- iV and (b) _E-_E-boson loop diagrams. 
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Fig. 2: Loop function -F(a;) of Eq. (5), vs x — m 2 N /M^y- 



need to have flavor changing neutral couplings (FCNC, 
absent in SM) to be relevant. But then they will face 
stringent constraints from /x — > ej and r — > £j. There- 
fore, if the 4th generation exists, Eqs. (l)-(4) will con- 
strain New Physics models. 

In the next section we first discuss the contributions 
involving neutral lepton N. In Sec. Ill we discuss the 
contributions involving charged lepton E. In Sec. IV we 
compare with the minimal supersymmetric SM (MSSM) . 
A summary is given in Sec. V. 



II. EFFECTS OF NEUTRAL LEPTON N 

The 4th generation neutral lepton N enters the one- 
loop diagram for muon g — 2 illustrated in Fig. 1(a). The 
boson can be the charged vector boson 14 7± , or charged 
scalar boson H ± , which we discuss separately. 

A. W + W~ N Loop Contribution 

With instead of N, this is the only contribution 
within SM. The contribution from a fourth generation 
lepton N has been considered before [15, 16]. We find 



a^W+W-N) 



472V 



F(x), 



(5) 



where 



rnj^/M^, is the lepton mixing matrix 



element, and 



du 



u(l - u)(2 - u)x + 2u 2 (u + 1) 
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We depict F(x) versus x in Fig. 2. We see that F(x), 
an Inami-Lim [17] loop function, is well-behaved and 
bounded, with F(l) = 17/12. However, this does not 
seem to be correctly rendered in Refs. [15] and [16]. In 
Ref. [15], though bounded, F{x) was not correctly eval- 
uated. Furthermore, the authors found strong enhance- 
ment near x ~ 1, which was used to put a bound on 



the fourth neutral lepton mass. However, from the inte- 
gral in Eq. (6) and the functional form of F(x), it should 
be clear that there is no enhancement near x — 1 (i.e. 
tojv ~ Mff). In Ref. [16], which is a study note for 
Ref. [18], the form of F(x) is again incorrect but still 
bounded. The author claimed that F{x) had a singu- 
larity at i = 1, and attributed this to the zero width 
approximation of the W boson propagator. Again, we 
see from Eq. (6), that there is no singularity for any x, 
and in any case, Tw should be irrelevant for such low 
scale processes. Our result therefore corrects some errors 
in the literature [19]. 

As we have already stated, if we replace by v^, 
we should recover the SM contribution. Using F(0) = 
5/3 in Eq. (5), together with V Uflll = 1, we get 

a SM(W+W-Vp) = 5G F ml/12V2n 2 . Furthermore, we 
find al M (^Z) = -(G F m 2 1 /6V2 7 r 2 )(l + 2sin 2 6W - 
4sin 4 #vy), where sin 2 9w = 0.23120. Combining the two 
together, and using 
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we get a^ M (l-loop Electroweak) = 195 x 10~ n , which is 
consistent with Ref. [20]. 

The SM exercise indicates that the 4th generation neu- 
tral lepton contribution has the right order of magnitude 
to contribute to Eq. (1). However, as seen from Fig. 2, 
the effect actually drops a bit from the massless result 
of SM as the mass of N becomes heavier. Furthermore, 
it is multiplied by the suppression factor | Vat^I 2 . From 
to 7v > 90 GeV [1], hence F(x) < 1.4, we see that \V N)i \ 
needs to be 0.7 or higher to reach within 2cr of Eq. (1). 
Considering the stringent constraint from Eq. (2), how- 
ever, this is clearly unrealistic. We conclude that the 
difference of Eq. (1) cannot come from the addition of a 
4th neutral lepton N. 



B. H + H~ N Loop Contribution 

It is unusual to consider both a 4th neutral lepton iV" 
together with charged Higgs H + . But since W + W~N 
contribution is insufficient for Eq. (1), we consider re- 
placing W + by the charged Higgs H + . This is the Two- 
Higgs-Doublet-Model (2HDM) with 4th generation lcp- 
tons. It is of interest to check whether one could gain 
from large tan/3 enhancement. 

For 2HDM-II (which occurs for MSSM), where up and 
down type quarks receive masses from different Higgs 
doublets, we find 
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Fig. 3: Loop functions f H +{ x ), 9h+( x )> 1h+{ x ) of Eq. (8), 
and h H +(x), qR+( x ) of Eq. (12) vsi = m 2 N /M^ 1+ . 
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functions in Eq. (8) are 
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/Mjj + and x M = m 2 l /M^ I+ . The loop 
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(1 — u) x + u 
x 2 log x ( 2 a; 2 + 5 a; — 1 
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We plot fji+(x), gH+( x ) an d Qh+( x ) i n Fig- 3. The 
q H +(x) term in Eq. (8) can actually be safely ig- 
nored, because one would need extremely large values 
of tan/3 to overcome the extremely small x^. However, 
we give a complete expression to check a previous re- 
sult in Ref. [21] for 3 generations. If we replace N be 
u^, one has f H +(0) = 0, g H +(0) = and q H +(0) = 
1/6, and Eq. (8) becomes a 2HDM - n (i/ + i?-^) = 

-(GFml/A^Xf, tan 2 /?/6. 

For the first term of Eq. (8), for 1 < x < 10 we have 
0-4 < fH+( x ) ^ 0.8, which is not particularly small. But 
because of the general |V/v/i| 2 suppression, an argument 
similar to the W + W~ N loop discussion suggest that this 
term can not give rise to Eq. (1). It is interesting that, 
because N has isospin +1/2, large cot/3 could lead to 
enhancement. If we take | Vn^ cot /3| 2 to be order 1 in 
the large cot/3 limit, and if mjy is large compared to 
, it could generate a finite contribution. But this 
contribution is negative, hence it is in the wrong direction 
for Aa M of Eq. (1). Furthermore, in the 2HDM-II (or 
MSSM), the UH°(h°) coupling relative to its SM value, 
m t /v, is given by cos a/ sin (3 (sin a/ sin (3). Large cot /3 
will make the coupling strength \g t iH° | ^ 1 or \gtih° I 3> 1 
and become nonperturbativc, which is not desirable. 



It has already been pointed, however, that there is a 
sufficient contribution in MSSM [20, 22, 23] coming from 
the large tan/3 region. We will give a brief assessment in 
Sec. IV. 

For 2HDM-I, where all quarks receive mass from the 
same Higgs doublet, we find 

,2HDM-Im+ 



\H + H~N) 



Gf m 2 ^ 
4V2V'^ 



; cot 2 f3[h H+ (x) - x M q H + (x)}, (12) 



with x and x^ as before 



h H +(x) 



da 



and 
u(l-u)(2- 



u)x 



x 2 (x 



(1 — u)x + u 
- 2) log x x(Jfx 2 



hx-h) 



is given in Eq. 
3. Analogous 



6(z - l) 3 



(13) 



(11). We plot h H +(x) 
to 2HDM-II, we have 
-H-Vp) = -{G F m 2 /A^/Tjx,, cot 2 /3/6, but 



while q H + (x) 
also in Fig. 

now everything is proportional to cot 2 (3. 

Similar to the 2HDM-II case, if we take \Vn^ cot /3| 2 to 
be order 1, it could generate a finite and positive contri- 
bution to Aa^. However, for 2HDM-I, the cot (3 enhanced 
Higgs couplings to tt are non-perturbative at large cot (3, 
which leads us to reject this possibility. 



III. EFFECTS OF CHARGED LEPTON E 



The 4th generation charged lepton E contributes to 
via i?-i?-boson loop diagrams, where the boson has 
to be neutral, and can be a scalar (h°), pseudo-scalar 
(A ), or vector (an extra Z'). But they need to possess 
flavor changing neutral couplings (FCNC). We discuss 
each case separately. 



A. EEh°, EE A Loop Contribution 

There is no /j,EH° coupling in SM. The same is true 
for the 2HDM-I and II, and fiEH , fiEh and f iEA° 
couplings are absent. This is because, by design [24], the 
charged leptons receive mass from just one doublet, and 
only one matrix needs to be diagonalized. In the so-called 
2HDM-III, this restriction is softened, and there exist two 
matrices rf^ and simultaneously for each lepton 

type. Note that in this model, by redefining <j>i, </>2 and 
r], £ simultaneously, which still leaves the Lagrangian in- 
variant, we may assume {(ff{) = v/y/2 and (</>!]) = with- 
out loss of generality, hence tan (3 is no longer a physical 
parameter. For a detailed analysis, we refer to Ref. [25]. 

To regulate the FCNC in face of stringent constraints, 
there is the ansatz suggested by Cheng and Sher [26] for 
the quark sector, i.e. all qiqjh / H° /A couplings have 
the same form 
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Fig. 4: Loop function F h o(x) of Eq. (15), vsi = m|/M ; 2 
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where Ay is of 0(1). This scheme can survive the rather 
critical constraint of K°-K° mixing, because y/mffml jv 
is extremely small. We extend it here to the charged 
lcpton sector with 4th generation, although it may not 
hold because the lepton mixing pattern seems different 
from those of the quarks. 

Note that CP-even Higgs bosons H°, h° give positive 
contributions to a M , but CP-odd A contributions are 
negative. Considering the positivity of Eq. (1), we may 
assume A is very heavy hence can be safely neglected. 
For sake of illustration, we set h° to be the lightest neu- 
tral Higgs, and assume no mixing between H° and h°. 
We then find 



Aa 



2HDM-III 



233 x 10" n F^o(a;), 



(15) 



where 



?/M 2 and we have taken A. y - = 1, and 



F h o(x) 



du ■ 



x log x x(x — 3) 



(.T-i) 3 2( x -iy 



(16) 



which is plotted in Fig. 4. 

There are other loops such as H + H~N to be con- 
sidered, but they are suppressed by m^/rriE and can 
be safely neglected. The suppression factor is ra^jwiE 
rather than (m^/ms) 2 because of the Cheng-Shcr cou- 
pling enhancement in Eq. (14). However, the LFV de- 
cay rates in Eqs. (2)-(4) give very stringent constraints, 
and need to be confronted. Here we use the formulas in 
Ref. [27]. Note that because a M and B{[i — > erf) come 
from loop diagrams of similar structure, their formulas 
are very closely related. After some organization, we have 



B 



2HDM-III 



(P -> erf) 



3a m e 
27r m u 



\F h °{x)f 



1.7 x 10~ 5 \F h a (x)\ 2 . (17) 



Let us consider first the case of r in the loop, which 
is the leading contribution with 3 generations, and was 



discussed in Ref. [28]. Eq. (17) becomes 

2J2HDM-III/ 



3a m e 
27r m„ 



,2 



S7)| 



3 gen. 



M 2 



log 



M 2 
M h° 



(18) 



Considering a factor of 2 uncertainty in Eq. (18), we still 
require M^o > 138 GeV in order to survive Eq. (2). We 
note that the MEG experiment can push the lower bound 
down to 530 GeV. 

Consider now 4 generations. Comparing Eqs. (1) and 
(15), we have Ff l o(x) = 0(1), and it seems that one 
could in principle bring about the difference with 4th 
generation under Cheng-Shcr ansatz. However, Eqs. (2) 
and (17) give F^o^x) < 10 -3 . From Fig. 4, we see that 
the <C M/jO is necessary, which is in conflict with data. 
Thus, the /it — ► &y constraint rules out the Cheng-Sher 
ansatz with 4th generation lepton in the loop. 

Even if we neglect the above shortfall, i.e. if we as- 
sume e decouples from Eh°, we still face the constraint 
on B(t — ► /17), i.e. 



B 



2HDM-III 



(r -»■ n"f) 



B(r^^ T )^^\F h o(x)\ 2 
3.6 x 10- 5 |F h o(a;)| 2 . 



(19) 



If we hold that m^/M/jO > 0.1 as reasonable, then t — > 
prj is again ruled out by Eq. (4) already. 

If the 4th generation is found, it seems that the Chcng- 
Sher ansatz can not hold for the lepton sector. 



B. EEZ' Loop Contribution 

For completeness, we discuss EEZ' contribution. Z' 
is the new gauge boson associated with an additional 
Abelian gauge symmetry U'(l) [29]. Because the typical 
constraint on the Z-Z' mixing angel 9 is < O(10 -3 ), 
we assume for simplicity that there is no mixing between 
Z' and Z. i.e. they are also the mass eigenstates. In 
terms of physical fields, the Lagrangian associated with 
the U'(l) gauge symmetry in the charged lepton sector 
is written as 



Cz>=gz>e a i > [e L ab L + e*R] e b Z'^ 



(20) 



where e^ b and e^ b arc the 4x4 chiral coupling matrices of 
Z' with charged leptons a, b = 1, 4 are flavor indices, 
and L(R) is the left(right)-handcd projection operator 
(1 =F 7 5 )/2. Because of the reality of the Lagrangian, 
e L (R) must be hermitian. 

After some calculation, we get the dominant contribu- 
tion to — 2 from EEZ 1 loop diagram, 



a,(EEZ') = & 



,1/2- 



'-Re(e^e§ M )/ z (X) 



c m (I 1 



ltE\ 



t nB\ 



Qz(x) , (21) 
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where x = m E /M^,, and x^ = rrvt/Mg,. The loop func- 
tions are 



fz{x) = 



qz{x) 



1 , 4u(l - u)x 1 / 2 + u 2 x z l 2 

du 

u x + (1 — u) 

3x 3 / 2 log x x 1/2 (x 2 + x + 4) 

(x-l) 3 + 2(x-l) 2 ' 

1 , u(l-u)(2-u) 

du — - 

q (1 — U) + U X 

x(2x — 1) log x 5x 2 + 5x — 4 



(x-1) 



4 



6(1- l) 3 



(22) 



(23) 



which are plotted in Fig. 5. 

As a check, we take the case of fifiZ. Then 
gz = e/ sin 9w coa9w, = 1/2 — sin 2 9w and e R = 
— sin 2 #w From Eqs. (22) and (23) we have fz{ x n) = 
2^fx~ and qz{x^) = 2/3. Hence, we find a^(ji[iZ) = 
-{G F m 2 /6^7r 2 )(l+2 sin 2 6^-4 sin 4 9 W ), which is con- 
sistent with Ref. [20], as mentioned in Sec. II-A. 

On the other hand, we can also get the contribution to 
the branching ratio of fi — > e"f from EEZ' loop diagram, 



2,ag%, 



Stt G\ m 2 Mi 



\f R e L I 2 



\f L f R I 2 
\ e eE e E f j,\ 



(24) 



For ttZ' in the loop, we find 



a^rrZ') = 



4tt 2 M 2 , 



which is given in Ref. [30], but with a sign error. A 
similar formula holds for E in the loop. It can be seen 
that if e^ E or e R E is zero (purely right-handed or purely 
left-handed), a^iEEZ') will become insignificant. This 
is because a spin flip is required. 

For sake of illustration, we take Mz> = ITeV, mg = 
250 GcV, gz> = 0.105 (predicted from a string model 
[31]), and denote e^ E = e R E = e^E and real, then 



a„(EEZ') 



620 x 10- n e 2 £ . 



(26) 



Furthermore, we assume e^ F = e R F = e eF and real, then 



2.47 x HT 2 e 2 B e 2 B . 



(27) 



Comparing Eq. (1) with Eq. (26), we need c^e = 0(1) 
for Z 1 to be the main contributor to muon g — 2. This 
is reasonable for a gauge interaction, as already stated. 
But one would need a model for why e^E is of order 
one. Furthermore, comparing Eqs. (2) and (27), we find 
e e E/£fj,E = O(10 -4 ) to satisfy /i — » e^f constraint. Hence, 
the model would not only need to account for e^E ~ 1, 
but e e E <C 1 as well. If this Z' is a main contributor to 
muon g — 2, it better not couple to electrons. 
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Fig. 5: Loop functions fz{x) and qz(x) of Eq. (21), vs x 
m%/M%,. 



For completeness, we discuss the contribution to the 
branching ratio of r — > ^7 from EEZ 1 loop diagram. 
Using a similar formula to Eq. (24), we find 



B(t -> [iv ll _v T \ 



3ag%, 



8nG 2 F m 2 M 2 , 



\fz(x)f 



\e R f L I 2 



\f L f R I 2 



(28) 



If wc take the same assumption and denote e^ E = ^e = 
e T E and real, then 

B z \r^^) ~ 1.52 x 10- 5 e 2 B e 2 £ . (29) 



(25) Comparing Eq. (4) with Eq. (29), taking e^s ~ 1 as 



before, wc find e rE ~ 10 is needed to satisfy r 



constraint. Thus, it does not seem likely that the Z' is 
the dominant source for muon g — 2. 



IV. COMPARISON WITH MSSM 

From previous discussions, we understand that the 
loop function (Inami-Lim functions) does not receive sig- 
nificant enhancement with heavy particle in the loop, 
so the coupling strengths become the crucial factors. 
For example, the magnitudes of electroweak contribu- 
tion a ll (W + W~ Vy) and a M (/i/iZ) both enter the inter- 
esting range 295(88) x 10~ n because they just involve 
gauge couplings, without any off-diagonal suppression. 
On the other hand, extreme smallncss of B(fi — ► ej) is 
predicted in the framework of SM, since there is no tree 
level FCNC, while loop effects are highly GIM suppressed 
by neutrino mass. Are other sizable electroweak contribu- 
tion possible? As we mentioned in Sec. II-B, there is 
such a mechanism [20, 22, 23] in the MSSM. 

Simply put, MSSM doubles the number of diagrams 
of SM. The corresponding loops to W + W~v^ and \x\iZ 
are chargino-chargino-j> M and /i-/i-neutralino respectively. 
Assuming mass degeneracy of supcrparticles, mHi ggs ino = 



TOWi] 



Mo 



Msxjsy, and in the large tan/3 limit (to 



5 



compensate the extra heaviness of Mrusy) [22], one can 
get a sufficient contribution. 

However, Ref. [20] used a different degeneracy condi- 
tion, TO c hargino = = MsuSYj which would inadver- 
tently send the chargino loop into a suppression region. 
Because it is a little subtle, we take a closer look. 

Following the formulas in Ref. [22] and under the con- 
dition TOchargino = = M S vsy, we have 



-9a:!/ 2 Re(C?*Cf )1 , (30) 



with x M = m£/Mf USY , Cf = {V^m^/v cos /3)(U x -) 2i 
and Cf = —g2(U x +)u, while U x + and U x - are unitary 
matrices and related by 



-Mr, 



g 2 v cos g 

Ga v/2 
sin /3 

7T~ 



fJ-H 



M SVSY U X+ Ul_. 



After some calculation, we get 



(31) 



(32) 



chargino-^ loop contribution is subdominant 



With Xf, = ml /M, 



Y ^ (10 7 ) typically, the chargino- 



V. SUMMARY 

In this paper, we consider the existence of 4th genera- 
tion leptons and discussed their impact on a M . In the SM, 
2HDM-I and II, the 4th generation seems irrelevant to 
the Aa^ puzzle because of the smallness of |Vjv>|. How- 
ever, this off-diagonal factor also protects these models 
from the stringent [i — ► ej and r — ► fi'y constraints. In 
the 2HDM-III, applying the Chcng-Sher ansatz with 4th 
generation to charged leptons, one has a strong conflict 
with B((J, — > ey) and even B(t — > ^7). Hence, if 4th gen- 
eration is found, the Cheng-Sher ansatz cannot hold for 
lepton sector. This may be reasonable since the lepton 
mixing pattern seems different from quarks. 

Our analysis illustrates why the well known SUSY 
mechanism is favored. Enhancement to a M and suppres- 
sion to B(fi — > &y) in the MSSM both bear similarities 
to the SM. It is interesting that the required large tan (3 
enhancement for the SUSY effect renders the negative 
contribution from H + H~N negligible, hence MSSM and 
4th generation can coexist. 
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